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In a recent work [Potter and Lee, Phys. Rev. Lett. 105, 227003 (2010)], it was demonstrated 
by means of numerical diagonalization that the Majorana end states can be localized at opposite 
ends of a sample of an ideal spinless p-wave superconductor with the strip geometry beyond the 
strict one- dimensional limit. Here we reexamine this issue, and study the topological quantum phase 
transition in the same system. We give the phase diagrams of the presence of Majorana end modes 
by using of Z2 topological index. It is found that the topological property of a strip geometry will 
change in an oscillatory way with respect of the sample width. 
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I. INTRODUCTION 

Recently, there is a growing interest in search- 
ing for Majorana fermions in the condensed matter 
physics community.^ Unlike traditional fermions, Ma- 
jorana fermions are their own antiparticles and are 
expected to obey non-Abelian statistics.—— Due to 
novel properties of the Majorana particles, it has 
been proposed that they can be used in realizing the 
fault-tolerant topological quantum computation pro- 
tected from local decoherence.— "— The early experi- 
mental efforts were focused on the 7 = 5/2 frac- 
tional quantum Hall state,— i^^ in which the Majorana 
fermions are expected to be charge e/4 quasiholes. 
Other proposed candidate materials include the p-wave 
superconductors ;2iitil^ in which the Majorana fermions 
are zero energy single-particle states trapped in vortex 
cores, and superfluids in the ^He-B phase.— »i^ Recent 
several new suggestions are proposed to look for Ma- 
jorana fermions, such as topological insulators proxi- 
mate to an s-wave superconductor, ^^•"•24 a semiconduct- 
ing thin film sandwiched between an s-wave superconduc- 
tor and a magnetic insulator^"— one-dimensional (ID) 
semiconductor-superconductor heterostructures based on 
quantum wires,— ^— noncentrosymmetric superconduc- 
tors with broken time-reversal symmetry,^^'^^^ and s- 
wave superfluid of ultracold atoms. '^^ Although exten- 
sive efforts have been made, direct evidence of Majorana 
fermions is still absent. Thus, it is still a challenging 
mission to produce and detect the mysterious Majorana 
particles. 

In his pioneering work, Kitaev^ found that Majorana 
particles can be localized at the ends of a ID p^^ -I- ipy su- 
perconducting wire. In a recent paper, '^^ Potter and Lee 
moved beyond the strict ID limit and explored further 
the notion of Majorana end states in quasi-lD p^ -\- ipy 
superconductors. They considered a spinless p-wave su- 
perconductor with finite width, and showed by means of 
numerical diagonalization that the Majorana end states 
can be localized at opposite ends of a sample with the 



strip geometry. We also note that a similar set of ideas on 
Majorana fermions in the quasi- ID wires was dicussed by 
Wimmer et ai2^ (to see Appendix C of the supplement of 
Ref. 35). Very recently, searching for Majorana fermions 
in multiband semiconducting nanowires attracts much 
attention and also opens a new path in the challenging 
field.—"— In this paper, motivated by the work by Pot- 
ter and Lee, '^4 we reexamine an ideal spinless quasi- ID 
Px + iPy superconductor, and study the topological quan- 
tum phase transition in the this system. Firstly, based on 
the tight-banding model, we give Majorana edge modes 
localized at boundaries of a sample with the sufficiently 
large width, and through analytical solutions in a model 
calculation for a strip of finite width, we find that two 
Majorana fermions on the two edges can couple together 
to produce a gap in the excitation spectrum under the 
periodic boundary condition along the longitudinal direc- 
tion. Then, we give the phase diagrams of the presence 
of Majorana end modes in quasi-lD sample by using of 
Z2 topological index. We find that the topological prop- 
erty of a strip geometry will change in an oscillatory way 
with respect of the sample width. 



The paper is organized as follows: in Sec. II, we shall 
give a square-lattice tight-binding model to describe a 
spinless p-wave superconductor. In Sec. Ill, we shall 
investigate the behaviors of Majorana edge states local- 
ized at boundaries of the sample with finite width, and 
the energy gap opened due to the coupling of two Ma- 
jorana fermions at opposite boundaries is given through 
analytical solutions in a model calculation^^ In Sec. IV, 
we shall calculate the Z2 topological index and plot the 
phase diagrams of the presence of Majorana end modes in 
the quasi-lD p-wave superconductors. By analyzing the 
phase diagrams, we shall then give the findings in the 
topological quantum phase transition in quasi- ID spin- 
less p-wave superconductors based on the tight-binding 
model. We shall summarize our conclusions in Sec. V. 



II. MODEL 

We consider a two-dimensional (2D) Kitaev model 
of spinless p-wave superconductors on a square-lattice, 
which is described by the following Hamiltonian^i^ 



H = 



n n 

j = l Q = l 

+tc]^^Cj+i^a + i |A| Cj^aCj + l^a + h.C. 



|A| Cj^aCj^a+l 



(1) 



where ri ^ creates an electron on site (j, a), t (> 0) is the 
hopping amplitude, /i is the chemical potential, | A| (> 0) 
is the p-wave pairing amplitude, and the lattice spacing 
is taken as unity. Here we assume a strip geometry in 
which the lattice site numbers are L along the x-axis 
direction and n along the y-axis direction (the sample 
width direction), thus which totals N — nL fermionic 
sites. First, one introduces a periodic boundary condition 



along the a;-axis direction, i.e., c'^^^ ^ 
the Fourier transform of the operator c- ^ 



-1,Q' 



and uses 



(2) 



where q = q^ is the wave vector along the x-axis, and 
— TT < q < TT. In terms of the new creation and annihila- 
tion operators cj, (q) and c^ (q) , the Hamiltonian ([T]) can 
be rewritten as 



H 



^ ^ { - (/i + 2t cos q)cl{q)ca{q) 

q a—1 
- [tC-i (q) Ca+l (q) + |A| Ca (q) Ca+1 {-q) 

+z |A| e^*«c„ (g) c„ (-g) + /i.e.] } . (3) 



Then, we define a set of the operators 720-1 (q) and 
72a (q) as 

72a-^i (q) = i [4 i-q) - Ca (<?)] , (4) 

l2a{q)^cl{-q)+Ca{q), (5) 

which satisfies the anticommutation relation 
{7^1 (9) ,7« {q')} = 2S.mr,5qg' and 7J, (q) = 7,„ (-g). In 
fact, 7m (0) is just a Majorana fermion operator due 
to 7m (0) = 7m (0)- I^ the basis of the news operators 
72a- 1 (<?) and 72a (q), the Hamiltonian (jS]) has the 
following form 



H = h'YYl ^^''^ (^) ^^ (~*) ^^ ("?) 



(6) 



where the elements of the 2n x 2n matrix B (q) are given 
as 




q/71 



FIG. 1: (color online). The excitation spectrum obtained 
by numerical diagonalization of the Hamiltonian with 
the parameters n = 100, A/t = 0.5. (a) fi/t = -2.0; (b) 
fi/t = 2.0. The (red) solid and dashed lines correspond to the 
edge modes, which are localized at opposite boundaries of the 
sample, respectively. 



-B2a,2a-1 = -S2Q-l,2a = -^ ^ 2< COS g, (8) 



-B2a,2a+1 — ^ -B2a+l,2a — —t — |A| 



(9) 



-B2a~l,2Q+2 — — -B2Q+2,2a-l — t — |A| , (10) 

and the else elements are zero. 

If we consider also the periodic boundary condition 
along the y-axis direction, the bulk excitation spectrum 
of a 2D square-lattice Kitaev model is given by 



E 



(q) = I [2t (cos qj; + cos qy) + fiY 

+A\A\^sm\,+sm\y)y . (11) 



1/2 



-B2a,2a = --620-1, 2a-l = -2i |A| siuq. 



(7) 



In general the spectrum (fTTj) has a gap, but the gap will 
close when fi = — 4i at q = {qx,qy) = (0, 0) or // = At at 
q = (±7r,±7r) or ^ = at q = (0,±7r) and q = (±7r, 0). 
Actually, |/i| = 4t is a phase transition line. For a strip 



geometry with the periodic boundary condition along the 
X-axis direction and the open boundary condition along 
the y-axis direction, if the sample width is sufficiently 
large (i.e., the lattice site numbers n is much larger than 
superconducting coherence length fo ^ V 1^1): ^ P^i'" of 
gapless chiral edge modes per boundary always present 
for —At < n < and < /i < 4t. In the former parame- 
ter region, the gapless points are present a.t q = q^ = 0, 
and in the latter at q = ±7r. The excitation spectra ob- 
tained by numerical diagonalization of the Hamiltonian 
^ with n — 100 are plotted in Fig. 1, where the chi- 
ral edge modes are shown inside the bulk gap and the 
energies of chiral edge modes are ±2 |A| sing. It is inter- 
esting to note that different from the continuum model 
the minimum of the edge dispersion shifts from q = 
to q = TT as one goes from negative to positive /i in 
the tight-binding model. This result can be understood 
based on the fact that if in the whole ID Brillouin zone 
(— TT < q < it) the gapless points are present at both q — 
and ±7r, the edge modes with the spectrum ±2 |A| sing 
would intersect Fermi energy (assumed inside the excita- 
tion gap) four times (an even multiple of two) which vio- 
lates the topological stability^i"— Hence the edge modes 
must merge into the excitation spectrum at finite q, and 
the gapless points occur either a.t q = or q = tt, but not 
at both.44 

As the sample width is decreased, due to the finite size 
effect, the wave functions of the edge modes on opposite 
edges overlap and mix, thus the spectrum reopens a gap. 
In Sec. Ill, we will focus on this issue. 



III. MAJORANA EDGE STATES IN THE 

BOUNDARIES OF p-WAVE 

SUPERCONDUCTORS 

Here we introduce a two-component field operator for 
each site a along the y-axis direction: 



da {q) 



l2a-l (q) 
72a (<?) 



(12) 



The Hamiltonian (jG]) can be expressed in the terms of two 
component creation and annihilation operators dj, (q) 
and da (q) as 



H = Y,{di (q) £da (q) + [di (q) fL+i (q) + h.c^ } 



q,a 



where 



£ = —2 |A| sinqaz — (/i + 2tcosq)ay, 



7"= -tay - i |A| cr^;, 



(13) 



(14) 



(15) 



with (7,; (i = X, y, z) being the Pauli matrices. We define 
the single particle states of the form 



'i'{q)=J2di{q)yJa\0}. 



(16) 



where |0) is the vacuum state annihilated by the opera- 
tor da ((?), and ipa = (i/'ai,'0a2) IS the two-component 
amplitude with index a. An open boundary condition 
Tpa^o = '4'a=n+i = is introduced along the sample 
width direction. The corresponding Schrodinger equa- 
tion is given by 

£ iq) ^a+f (q) iba+l + P (q) V-a-l = £ (?) i^a- (17) 

Note that Eq. ([T7| is analogue to the equation given in 
studying edge modes of topological insulators based on 
the tight-binding model. ^'^"''^ 

We study solutions of the form tpa ^ A"?/) with A be- 
ing a general complex number 4^ Firstly substituting the 
ansatz into Eq. (fT7|) . we have 



'£iq) + Xfiq) + X-'fHq)\i' = eiq)^J 



(18) 



Thus the secular equation gives an equation about A and 

eiq) 



{t' - \Af) 



't^ - I AH z^-f 2t(^ + 2icosg)z 
= e^ -A\l^\^ {I + sm^q)-{^i + 2t cos qf , (19) 

where z = (A + A^^). Note that, if A is a solution, so 
is A~^. Thus for every exponentially increasing solution, 
there exists another one which exponentially decreases4i 
Four roots of A are given by 



Ai,: 



z± 



zl-A 



X 



z± 



3,4 



(20) 



where z± are two roots of z in Eq. (|T9|). We define 
two of four roots of A as A± which satisfy the relation 
|A±| < 1, and then the other two roots of A are A^^ 
satisfying |A^^| > 1. 

Now we consider the edge states. For convenience, a 
symmetric boundary condition is used in the following, 
that is V'a=-F/2 = "05=^/2 = 0. Here, y = n+ 1, and the 
new index a = a— {n + 1)/2. Thus for the wave function 
"05 — ("051, "0152) we have an analytical expression' 



48 



ip&i = C+/+ [q, a) + c_/_ (g, a) , 



ipa2 = d+f+ (q, a) + d-f- (g, a) 



(21) 



(22) 



where 



cosh(A+a) cosh(A_Q;) 
f+{q,a)^ 7^ ^ 7^ ^, (23) 



cosh 



(~X+Y/2] cosh (A_r/2] 



sinh I X+a 



f~ (q, a) = 



sinh f A_(5 j 



sinh(A+r/2) sinh(A_y/2 



(24) 



with A± = lnA±. The nontrivial solution for the coefS- 
cients c± and d± in the wave functions leads to a secular 
equation 

|2 



t-" - lAH f cosh^ A+ + cosh^ A_ 



2|A| 



2t^ coshA+ coshA- - T |A| sinhA+ sinh A_,(25) 



with 



T 



tanh(A+y/2) tanh(A_y/2 



tanh 



(~X-Y/2\ tanh (~X+Y/2 



(26) 



Equations ([20| and (|25|) give the energy dispersions and 
the values of two characteristic quantities A± — e^± . 

When the sample width is sufficiently large, i.e., in the 
limit of n — >■ oo, we can find two chiral edge modes with 
energy 



£±((?)=±2|A|sin(7, 



(27) 



and 



A± = 



-(^ + 2icosg)± J{fi + 2tcosqf - W^ +4|A|^ 



(28) 

For nonzero q, the eigenstates for e± > become concen- 
trated on one edge or the other, depending on the sign of 
g, that is the left-moving and right-moving edge modes 
are localized at the opposite boundaries, respectivelyi^ 
Equation ([27l) is consistent with the result by numerical 
diagonalization (to see Fig. 1). Interesting to note the 
cases oi q — (for ^ < ) and tt (for ^ > ), one has 
zero energy mode e {q = O/tt) = which corresponds to 
the Majorana edge state in p-wave superconductors.^ 

Actually, there are no exact e = modes for a finite 
width sample. One of the key features for the solution 
of a finite width is the gap Eg opening for the energy 
dispersion of the edge state. From Eqs. (PU]) and (PSJ). 
one has 

e'^ -ilAfsiii^q 

= 2t{ii + 2t cos q) (cosh A+ + cosh A_ ) + (/i + 2t cos q) 

-h4i2 cosh A+ cosh A_ - 2r |A|^ sinh A+ sinh A_, (29) 

thus it can be found that a finite energy gap Eg at q — 
(for the case of /i < ) is approximately 



E„ 



\ 



2fi in + 4t) 



f2-|Ar 

and at g = TT (for the case of ^ > ) 



-n/io 



E„ 



\ 



2n {fi - At) 



t^ 



i/lo 



(30) 



(31) 



where ^g"^ = min(|ln|A+|| , |ln|A_||), and lo denotes the 
localization length of the edge mode. The gap Eg will 
decrease exponentially with increase of the width of the 
strip. 



IV. MAJORANA END STATES IN QUASI-ID 
p-WAVE SUPERCONDUCTORS 

In a recent work by Potter and Lee;^ two well-isolated 
Majorana end states locahzed at opposite ends of the 
strip geometry have been obtained by numerical diago- 
nalization of the Hamiltonian ([T]). Here we will give the 
phase diagrams of the presence of Majorana end modes 
in quasi- ID p-wave superconductors by using topological 
arguments due to Kitaev,^ To this aim, we consider the 
2n X 2n matrix B (q) in the Hamiltonian ([B|) . The matrix 
B is an antisymmetric matrix when q is equal to zero 
or TT, such that we can calculate the PfafRans Pfi? (0) 
and PiB (tt). The topological property of the system de- 
scribed by the Hamiltonian ^ is characterized by a Z2 
topological index (Majorana number) M: 



M ^ sgn [PiB (0)] sgn [PiB (tt)] ^ ±1, 



(32) 



where -1-1 corresponds to topologically trivial states and 
— 1 to topologically nontrivial states (i.e., the existence 
of zero mode Majorana end states) i^ '^^i'^^'^^ 

For the simplest case, there is only one lattice site along 
the y-axis direction (i.e., n = 1). This case is just the 
Kitaev original model.- Two 2x2 antisymmetric matrices 
are 



Bn^i (O/tt) 



^l±2t 

in ± 2t) 



(33) 



and PfB„=i (O/tt) = iJ,±2t, where "-h" and "-" corre- 
spond to the cases oi q = and tt, respectively. The 
Majorana number for the case of the strict ID limit is 
given 

7W„=i=sgn(/i + 2t)sgn(/i-2f), (34) 

thus we have the topologically nontrivial condition 

2|t|>H (|A|^0). (35) 

The above Eq. ((35)) is just the result by Kitaev,-^ who 
demonstrated for a long open chain (in the limit of 
L — >■ 00) there are zero energy Majorana end states local- 
ized near per boundary point under the condition ([35t . 
However, if the chain length L is finite, there is a weak 
interaction between two unpaired Majorana fermionsj^ 
Following the method given in Sec. HI, we can also ob- 
tain the energy gap opened due to the finite size effect 
approximately as 



A.= |A| 



\ 



2 (4^2 - f^2) 



~L/Iq 



i2-|A|^ 

where Iq = min(|ln |a;+|| , |ln|a:_||), with 



(36) 



x± 



2(i+|A|) 



(37) 



Here, Iq indicates the localization length of the Majo- 
rana end states. Thus, the energy gap Ag vanishes as 

exp ( —L/lo 1 in an open chain.— 

For the case of n = 2, the lattice site numbers along 
the y-axis direction are two. Two 4x4 antisymmetric 
matrices are 



Bn=2 (O/tt) 

n±2t i-|A|' 

-(/i±2t) -{t + \A\) 

i+|A| fi±2t 

■(t-|A|) -(^±2<) . 



.(38) 



The direct calculation yields the Pfaffians Pfi3„^2 (O/tt) 



PiBn=2 (0/^) = {fi± 2ty + |Ar - f 



(39) 



For the lager lattice site numbers n (> 3), Pfi3„ (O/tt) 
can be also calculated analytically, and we obtain a re- 
cursion relation 

PfB„ (0/^) ^ a±PiBn-i (O/vr) + bPiBn-2 (O/tt) , (40) 

where a± = iJ,zL2t and 6 = |A| — i^. We further solve 
Eq. (|40)) . and give an analytic formula for PiB^ (O/tt) 



PfB„ (0/^) = 



ir^'-^') 



46 



(41) 



where 



n 



a± 



46 



a± 



r2 



46 



(42) 



According to the Pfaffians Pfi?„ (O/tt) one can compute 
A^ as a function of the physical parameters, and then 
plot the phase diagram showing a sequence of topolog- 
ical phase transition for different the lattice site num- 
bers n. Figures 2 and 3 plot the phase diagrams for 
the even and odd lattice site numbers n along the y- 
axis direction, respectively. The phase diagrams of this 
tight-binding model have the symmetry on positive and 
negative /x values, thus here we only plot on negative /i 
values because the other part on positive ji values is a 
mirror image. However, this /i — > — /i symmetry is not 
generic to models with say, next-nearest-neighbor hop- 
ping or next-nearest-neighbor pairing. 

We now analyze these phase diagrams. Firstly, it is 
observed that for the given value of |A| /t the topolog- 
ically nontrivial and trivial phases alternate with the 
variation of the value of ^/t. As an example, the case 
of |A| /t = 0.1 is shown in the Figs. 2 and 3. From 
the Figs. 2 and 3, it is shown that the phase diagrams 
have different properties for the even and odd lattice site 
numbers n. Such the odd-even effect should be an arti- 
fact of the tight-binding model and may not be generic 
to all quasi-lD p-wave supercondurtors, since the order 
parameter for Cooper pairing are presented as a model 
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FIG. 2: (color online). Phase diagram for the quasi-lD p-wave 
superconductor model as a function of the p-wave pairing am- 
phtude and chemical potential for even lattice site numbers 
n along the y-axis direction. "N" denotes the topologically 
nontrivial region in the present of zero mode Majorana end 
states and "T" denotes the topologically trivial region with- 
out zero mode states. When |A| /t — 0.1, the solid (red) lines 
and dotted (blue) lines guide the values of /i/t corresponding 
to the topologically nontrivial and trivial phases, respectively, 
(a) n = 2; (b) n = 4; (c) n = 6; (d) n = 8; (e) n = 10. 



parameter, not given by a self-consistent mean field cal- 
culation. Even in a realistic material, the oscillation of 
the critical temperature with the monolayer number of 
the sample was observed experimentally as quantum con- 
finement effect.— It is not clear what the real meaning 
of odd and even n in this model calculation, which may 
have no counterpart in the continuous model as a long 
wave length limit. However, the oscillation of the pres- 
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FIG. 4: (color online). Excitation spectrum Eraj |A| for a 
square-lattice Kitaev model obtained by numerical diagonal- 
ization of the Hamiltonian ([ij with parameters: L — 100, n — 
10, I A| /t = 0.1. m labels eigenvalues of the Hamiltonian ((TJ. 
Here only the partial low excitation energies are presented, 
(a) /i/t = -1.0 (£'i/|A| = 1.09 X 10"^); (b) fi/t = -1.5 
{El/ |A| = 0.1116); (c) ^I/t = -2.0 {Ei/ |A| = 1.54 x 10"^); 
(e) fi/t = -2.5 [El/ |A| = 0.1595). 
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FIG. 3: (color online). Phase diagram for the quasi- ID p-wave 
superconductor model as a function of the p-wave pairing am- 
plitude and chemical potential for odd lattice site numbers n 
along the y-axis direction. "N" denotes the topologically non- 
trivial region in the present of zero mode Majorana end states 
and "T" denotes the topologically trivial region without zero 
mode states. When |A| /t — 0.1, the solid (red) lines and dot- 
ted (blue) lines guide the values of /i/f corresponding to the 
topologically nontrivial and trivial phases, respectively, (a) 
n = 3; (b) n = 5; (c) n = 7; (d) n = 9; (e) n = 11. 



ence and absence of the zero mode still clearly exhibits 
in the continuous model as quantum confinement effect. 

According to the topological argument above, the 
topologically nontrivial phase corresponds to the exis- 
tence of the zero energy mode, which can be justified by 
the excitation spectrum obtained by the numerical di- 
agonalization of the Hamiltonian (jlj . According to the 



FIG. 5: (color online). Excitation spectrum Em/ |A| for a 
square-lattice Kitaev model obtained by numerical diagonal- 
ization of the Hamiltonian ([T} with parameters: L = 100, 
|A[ /t = 0.1, n/t = —1.1. 771 labels eigenvalues of the Hamil- 
tonian IT]). Here only the partial low excitation energies 
are presented, (a) n = 6 (Si/ |A| = 0.0966); (b) n = 7 
{El/ |A| = 8.21 X lO"*^); (c) n = 8 {Ei/ \A\ = 7.39 x lO"''); 
(e) n = 9 (£i/|A| =0.1102). 



work by Potter and Lee,'^^ for sufficiently long samples 
(L ^ e"/^"), one can find spatially isolated zero energy 
Majorana end states localized at opposite ends of the 
sample with the strip geometry. To see Fig. 2(e) (with 
n = 10), we can find that if |A| /t is fixed at 0.1, the cases 
of fi/t = —1.0 and —2.0 correspond to the topologically 
nontrivial phase; the cases oi fj,/t = — 1.5 and —2.5 to the 
topologically trivial phase. For the above four parameter 
cases, the excitation spectra of the Hamiltonian ([IJ with 
L = 100 and n = 10 are plotted in Fig. 4. Numerical 



diagonalizations show that the zero energy modes indeed 
exist in the cases of ^x/t = — 1.0 and —2.0; there are no 
zero energy mode in the cases of /i/f = —1.5 and —2.5. 
The conclusions are consistent with the resuhs by Potter 
and Leej2i 

We note also an interesting property from the phase 
diagrams. For some given parameters, the topological 
property of a strip geometry will change in an oscillatory 
way with respect of the sample width. For instance, in 
Fig. 2, when the parameters |A| /t = 0.1 and ji/t = 
—2.0 are fixed, the topologically nontrivial and trivial 
phases alternate when the lattice site numbers n (even 
numbers) along the y-axis direction change from n = 2 
to 71 = 10. While in the cases of the odd n, ^x/t = 
—2.0 is a topological phase transition line (to see Fig. 3.) 
As an other example, we take the parameters |A| /t = 
0.1, /i/t = —1.1, and then from the Figs. 2 and 3 it is 
presented that the topologically nontrivial phases exist 
in the cases of n = 2, 4, 5, 7, 8, 10, and 11; while the 
cases of 71 = 3, 6, and 9 correspond to the topologically 
trivial phases. In Fig. 5 the excitation spectra for four 
different sample widths (ii = 6, 7, 8, and 9) are shown. 
For the case n = 7 and 8, there are zero energy modes 
in the given parameters; while the zero energy modes 
disappear in the case n = 6 and 9. Thus the topological 
property obtained by the phase diagrams are consistent 
with the results by numerical diagonalization. 



V. CONCLUSIONS 

In this paper, we investigate crossover from Majorana 
edge to end states in an ideal spinless quasi-lD p-wave 



superconductor based on the tight-banding model. We 
found the existence of Majorana edge modes when the 
sample width is the sufficiently large, and then through 
analytical solutions in a model calculation two Majorana 
fermions at the two edges can couple together to pro- 
duce a gap in the excitation spectrum for a strip of finite 
width under the periodic boundary condition along the 
longitudinal direction. We calculate Z2 topological index 
in the quasi-lD sample, and plot the phase diagrams of 
the presence of Majorana end modes. By analyzing the 
phase diagrams, we find that for some given parameters 
the topological property of a strip geometry will change 
in an oscillatory way with respect of the sample width. 
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